Spin Squeezing under Non-Markovian Channels by Hierarchy Equation Method 
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We study spin squeezing under non-Markovian channels, and consider an ensemble of N inde- 
pendent spin-1/2 particles with exchange symmetry. Each spin interacts with its own bath, and 
the baths are independent and identical. For this kind of open system, the spin squeezing under 
decoherence can be investigated from the dynamics of the local expectations, and the multi-qubit 
dynamics can be reduced into the two-qubit one. The reduced dynamics is obtained by the hierarchy 
equation method, which is a exact without rotating-wave and Born-Markov approximation. The 
numerical results show that the spin squeezing displays multiple sudden vanishing and revival with 
lower bath temperature, and it can also vanish asymptotically. 

PACS numbers: 03.67. Mn, 03.65. UD, 03.65. Yz 



I. INTRODUCTION 

Spin squeezing has attracted much attention for 
decades An important application of spin squeez- 

ing is to detect quantum entanglement |Ql— Till] . As a 
multipartite entanglement witness, spin squeezing is rel- 
atively easy to be generated and measured @, E2l4l4j ]. 
Many efforts have been devoted to find relations be- 
tween spin squeezing and entanglement [lrE], Il5l - fl7j . An- 
other application of spin squeezing is to improve the pre- 
cision of measurements. For example, spin squeezing 
plays an important role in making more precise atomic 
clock B, IB, EH EH and gravitational- wave interferome- 
ters [2fj|422| . and so on. 

Spin-squeezed states are useful resources for quantum 
information processing. However, in practice, decoher- 
ence is inevitable and harmful to spin squeezing and en- 
tanglement [23:- 29]. Generally, when the system-bath 
coupling strength is weak enough, the decoherence is 
studied by using the master equation method, which is 
derived by employing the Born approximation (23l . l24j . 
Besides, the Markov approximation can be applied if the 
time scale of the bath is much shorter than that of the 
system. To overcome the above approximations, a set of 
hierarchical equations were established by Tanimura et 
al (30l - [3r3 |. It provides an exact way to obtain the re- 
duced dynamics of system (37l.[38|. However, for numeri- 
cal reasons, it is hard to treat systems with large number 
of particles straightforwardly. Here, we show that for 
the open system we consider, we can reduce the multi- 
particle dynamics to the two-particle one, and then we 
efficiently use the hierarchy equation method to make 
numerical calculations. 

As we know, spin squeezing is a multipartite entangle- 
ment witness. Reference [39| has shown that for a many- 
particle system with exchange symmetry, the spin squeez- 
ing parameters of the total system can be expressed in 
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terms of local expectations and correlations. Here, we 
consider such an ensemble of N independent spin-1/2 
particles. Each particle interacts with its own bath, and 
the baths are independent and identical. Thus, the ex- 
change symmetry is not affected by the decoherence, and 
the spin squeezing parameters of the open system can 
be expressed by dynamics of the local expectations and 
correlations. For the system under consideration, we find 
that the dynamics of any two particles is governed only 
by the local Hamiltonian of the two particles and their 
baths. Then, we use the hierarchy equation method to 
calculate the dynamics the the local expectations and 
correlations. Reference ;39j has also shown that the spin 
squeezing has close relation with pairwise entanglement if 
the state of the collective spin system lies in the J = N/2 
sector, where J is the collective angular momentum of the 
system. Therefore, since the state of the system will not 
lie in J = N/2 sector anymore under decoherence, the 
ability of spin squeezing in detecting pairwise entangle- 
ment needs to be further studied and clarified. 

This paper is organized as follows. In Sec. II, we in- 
troduce the Hamiltonian and the initial state of the open 
system. The definition of the spin squeezing parameters 
is given in Sec. Ill, and we also discuss the symmetry 
of the system and reduce the multi-qubit dynamics into 
the two-qubit one. In Sec. IV, we introduce the hierarchy 
method and give a alternative form of the hierarchy equa- 
tion. We numerically calculate spin squeezing parameters 
and the rescaled concurrence of the open system under 
decoherence and compare their behaviors in Sec. IV. At 
last, a summary is given in Sec. V. 



II. HAMILTONIAN AND INITIAL STATE 

The system we consider is an ensemble of N inde- 
pendent spin-1/2 particles with exchange symmetry, and 
each particle interacts with its own bosonic bath. The N 
baths are independent and identical. The Hamiltonian 
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of the total system is (ft = 1) 

H = Hs + Hb+Hsb 



N 



E w \ 



a=l 
JV 



^ ^ Qak&ax 
a — 1 k 



(4 + «fc) 



(i) 



where the first term is the Hamiltonian of the system 
with aka(a = x, y, z) the Pauli matrices for the fc-th spin 
and luq the frequency for all qubits. The second term 
describes the bosonic bath, where bk and b\. are the cre- 
ation and annihilation operators of the fc-th mode with 
frequency ujk ■ The system-bath coupling is characterized 
by the third term with g a k the coupling strength for qubit 
a. Here, we study N independent baths, i.e., the bath 
can be divided into N parts and g a k is only non-zero 
when mode k belongs to the a-th part. 

The initial state of the total system is set to be a prod- 
uct state 



Pt(0) =ps{0)®Pb{0), 



(2) 



where ps(0) is a spin-squeezed state and pb(0) is a ther- 
mal state given by 



pb(o) = n 



cxp(-/3w fe c4a/ c ) 



Z k 



(3) 



with the inverse temperature j3 — l/(fcsT) and partition 
function = Trexp(— j3uka\a.k) for mode k. In this 
paper we take kg = 1- 

We choose the initial state as a standard one-axis 
twisted state [l[ 



l*(0)> 



1 1- 1) 



with 



1 N 

J a = -Z y^O-fcc 

k=l 



(4) 



(5) 



the total angular momentum operators. This state is 
prepared by the one-axis twisted Hamiltonian H = X^xi 
with the coupling constant \ , an d $ = the twist 
angle. For our case, the system of N spin-1/2 behaves 
like an effective large spin N/2. 



III. SPIN SQUEEZING AND REDUCING THE 
MULTI-QUBIT DYNAMICS INTO A 
TWO-QUBIT ONE 

In this section, we give the definitions of two spin 
squeezing parameters. By discussing the symmetry of 
the open system under consideration, we know that the 
spin squeezing can be expressed by the local expectations 



and correlations. Since we can reduce hte multi-qubit dy- 
namics into a two-qubit one, the spin squeezing can then 
calculated by the dynamics of the local expectations and 
correlations. 



A. Spin squeezing definitions 

There are various measures of spin squeezing related 
to various inequality criteria [U-III, S HI j and we consider 
two of them as follows: 
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N 



(J 2 ) 



N 
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(6) 
(7) 



Here, the minimization in the first equation is over all 
the directions denoted by _L, which are perpendicular to 
the mean spin direction (J)/\(J)\. A m i n in the second 
equation is the minimal eigenvalue of the matrix 



where 



Iki = Cfcz - (Jk)(Ji), k,le{x,y,z} 



is the covariance matrix and 



Cfej — -^(JiJk + JkJl) 



(8) 



(9) 



(10) 



is the global correlation matrix. The parameters ^jLj and 
£t were defined by Kitagawa and Ueda [lj, and Toth et 
al. [H, respectively. If ^ < 1, spin squeezing occurs, and 
we can safely say that the multipartite state is entan- 
gled@,i. 

From the definitions, we know that the spin squeezing 
parameters are based on the expectations and correla- 
tions of the collective operators. For the limitation of 
the hierarchy equation method, it is hard to calculate 
the decoherence of the many-particle system straightfor- 
wardly 



B. Simplification of the spin squeezing parameters 

Since the baths are independent and identical, the ex- 
change symmetry is not affected by decoherence. There- 
fore, the global expectations or correlations of collective 
operators can be written as [39j 



TV 



(Jl) 
[Ja, Jp} + ) 



N N(N — 1) 
T + 4 
N(N - 1) 



\&laO'2a), 



(11) 

(12) 



[via, 0-2/3]+), (a ^ /3),(13) 
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which only depend on the expectation values of the local 
Pauli operators, e.g., (ai a a 2 ^) an d (<Ji a ). 

The initial one-axis twisted state we use here has a 
parity symmetry leading to (J x ) = (J y ) = 0, namely the 
mean-spin direction is along the z-axis. Moreover, the 
mean-spin direction do not change during decoherence. 
The proof is given as follows. 

The Hamiltonian (flj displays only one symmetry, i.e., 
the parity symmetry. The parity operator is given by 



n = ni®n 2 



(14) 



where Af = J z + N/2 describes the numbers of excitations 
of up spins. Obviously, we have 



T1HII = H, 

nip s (o)n! = Ps (o), 

n 2 p B {o)n 2 = PB (o), 

n PT (o)n = p T (o), 



(15) 
(16) 
(17) 
(18) 



namely, the Hamiltonian and the initial state have a fixed 
parity. Since the exchange symmetry leads to (J x ) = 
N(ai x )/2, we obtain 

(a lx ) = Tr[a lx U(t)p T (0)tf(t)] 

= Tr{<7i x n[nc7(t)n] [n PT (o)n] [nzj+(t)n]n} 

= Tr[a lx nC/(t)p T (0)C/t( t )n] 
= Tr[lla lx TlU(t)p T (0)uHt)] 

= -{*!*), (19) 

which leads to (J x ) = 0. Similarly, (J y ) — (J y J z ) = 
(JxJz) = can be proved. Therefore, during the evolu- 
tion the mean spin direction is always along the z-axis. In 
this case, the spin squeezing parameters reduce to 0, [H[ 

&u - l + 2(iV-l)«a 1+ a 2 _>-|<<7 1 _a 2 _)|), (20) 



min {^c/^ 2 } 



(l-l/N)(ai-a 2 ) + 1/N' 



where 



? 2 = 1 + (N - 1) ((a lz a 2z ) - (a lz )(a 2z )) 



(21) 



(22) 



For convenience, hereafter we use 

C2=max(0,l-$), fce{KU,T} (23) 

to characterize spin squeezing. With the above definition, 
spin squeezing occurs when ^ > 0. 

Now we only need to calculate the dynamics of the lo- 
cal expectations and correlations of the spins, the spin 
squeezing parameters are greatly simplified. However, 
we still have to use the dynamics of the density matrix 
of the system to calculate the local expectations and cor- 
relations. 



C. Reducing the multi-qubit dynamics into a 
two-qubit one 

Now we prove that we can reduce the multi-qubit dy- 
namics into a two-qubit one. Generally, we consider a 
system written as follows 



N 



SB- 



(24) 



i=l 

(i) (i) 

Hg is the Hamiltonian of one particle, H B is the bath 

(i) 

Hamiltonian, and the couplings are expressed by H SB . 
Obviously, each of the particles interacts with its own 
bath. The particles do not have interaction with each 
other, and the baths are independent. Equation ([1} be- 
longs to this case. 

The time-evolution operator of the total system can be 
written as 

U(t) = e- lHt = H e- iH ^ = J] Ui (t), (25) 



where Ui(t) — e ** . Then, the total density matrix at 
time t is given by 



PT (t) = U(t)p T (0)U^t), 
which can be formally written as 

p T (t) = U(t)pr(P)U*(t) 



(26) 



(27) 



Here we assume that the initial state is a product state 
written as 



Pt(0) =Ps(0)®Pb(0). 



(28) 



By tracing out the baths and N — 2 particles of the sys- 
tem, we obtain the reduced density matrix of any two 
particles 



pf(t) 

Tl {Bi,2} 
Tl "{Bi,2} 



N 



N 



= Tr 



{Bl,2> 



\i=l i=l / 

/ 2 2 N 

^{s 3 ... n b 3 ... n} n u *(*)^(°)iNw 

\i=l i=l / 

2 2 



(29) 



where the second equality follows from the fact (see Ap- 
pendix A) 

Tr 2 [{A 1 ®A 2 )p 12 {B 1 ®B 2 )\ 
= Tr 2 [At ® (B 2 A 2 )p 12 (Bi ® I 2 )} , (30) 
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and the last equality is obtained by substituting the ini- 
tial product state ([2g]). pf (0) = Tr {S3 N }Ps(0) and 
Pb (0) = Tr{ Ba jyjPsCO) in the equation are the reduced 
density matrices of the initial state for the system and 
baths respectively. 

The above equation (|30| tells that the evolution of any 
two particles is governed only by the local Hamiltonian 
of the two particles and their baths. It is noted that we 
can reach this conclusion even when the initial state of 
the system or the baths are entangled states. Therefore, 
the multi-qubit dynamics reduces to the two-qubit one. 
Then we use the hierarchy equation method to calculate 
the two-qubit reduced density matrix of the system, and 
the dynamics of the local expectations and correlations 
in Eqs. (|20 |) - ([22l) can also be obtained. 

Here we emphasize that we obtain this conclusion with- 
out using exchange symmetry, which means that the par- 
ticles are not necessarily identical, and so do the baths. 
Also, the proof can be easily extended to any finite num- 
ber of particles. 



IV. HIERARCHY EQUATIONS AND INITIAL 
TWO-QUBIT REDUCED DENSITY MATRIX 

To start with the numerical calculations, we introduce 
the hierarchy equation method [H, [37[ and discuss the 
spin squeezing parameters of the initials state in this sec- 
tion. For comparison, the definition of a rescaled concur- 
rence is also given. 



A. Hierarchy equations 

We choose the Drude-Lorentz spectrum, 

2 WA7 



J(u>) 



9 i 1' 



where 7 represents the width of the spectral distribution 
of the bath mode and A can be viewed as the system-bath 
coupling strength. The bath correlation function for the 
bath operator 



k 

is given by [37{ 

00 

(B a (t)B a ( T )> = Y,Cne- Unlt 

where 



iui k t 



n=0 



Ink, 

is the A:-th Matsubara frequency, and 



Ck = —5 — 5 ^(l-^feo) + A7 



cotl T 



(32) 



(33) 



(34) 



<5fco (35) 



are the expansion coefficients. 

With the Drude-Lorentz spectrum, the hierarchy equa- 
tions becomes 

PH = - [iHg + (rti + n 2 ) ■ v] pa 

2 M 

-i X! X! nQfc ( c kVaPn-e ak - c*kPn-s ak V a ) 



a=l k=0 
2 M 

a=l k=0 



(36) 



where 



n= (ni,n 2 ) = (mo, uim, «20, «2m) (37) 

is a 2(M+ l)-dimensional vector, a concatenation of two 
(RI + l)-dimensional vectors n\ and ni. The vectors v = 
(yo, ...I'm) and e a k are defined as 2(M + l)-dimensional 
vectors with only 1 in the ak place and 0s in other places. 
Note that this equation is slightly different and essentially 
the same as that given in Ref. [37j . 



B. Initial two-qubit reduced density matrix 

To solve Eq. (|3"6"|). we need the initial state. Since the 
mean spin of the initial state Q is along the z-direction, 
the two-qubit reduced density matrix can be written as 
a block-diagonal form Q , 



m = 1 v+ u ® r V ) . ( 38 ) 



(31) in the basis {|00), |11), |01), |10)}, where 



v ± = (1 ± 2(a u ) + (<ru<T2z)) /4, (39) 

w = (l-{a u a 2z })/4, (40) 

u = (ai-o-2-), (41) 

y = (a 1+ a 2 -). (42) 

We notice that if (ai + a 2 -), (<ti_(T2-), (ci z ), and 
(<7i z (T2z) are known, the density matrix is determined. 
For the one-axis twisted state, we have 0] 



N-l 



(a,) = -cos 



((Ti_(T2_) 



1 — cos 



2j ' 

N - 2 e) 

N—2 a\ 



[ 1 — COS 



N-2 , 



--Bin I - I cos 



N-2 



(43) 
(44) 
(45) 

(46) 
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Employing the equations above, we obtain the initial two- 
qubit reduced density matrix in Eq. (|38[) . Then we use 
Eq. (|36| to calculate the dynamics of the reduced density 
matrix numercially. 

Meanwhile, we can also use the Eqs. p3)) - (l46l) to dis- 
cuss the spin squeezing parameters for the initial state. 
For the initial state we obtain 



Cku(O) = Ct(0) = 



4 

2JV-4 



(1 — COS 

1/2 



JV - 2 e) 2 + 16sin 2 



1 + COS 



N-2 , 



, (47) 



which implies that the two spin squeezing parameters for 
the initial state coincide. 

It is known that the spin squeezing has close relation 
with concurrence if the state of the collective spin system 
lies in the J = N/2 sector [3^|, such as the initial state 
of the system. During the decoherence, the state of the 
system does not lie in J — N/2 sector anymore. It is 
necessary to compare the behaviors of spin squeezing and 
pairwise entanglement. 

The concurrence is defined as 40] 



C = max(0, Ai — A 2 — A 3 — A4), 



(48) 



where Ai > A2 > A3 > A4 arc the square roots of eigen- 
values of pp. Here p is the reduced density matrix of the 
system, and 



P= {<7y®<7y)p*{Vy®<7 y )i 



(49) 



where p* is the conjugate of p. 

For the reduced density matrix of (|38[) , the concurrence 
is given by [Il[ 



C = 2max{0, \u\ 



y,y 



\/V+V-} 



(50) 



Therefore, we can also obtain the concurrence of the ini- 
tial state by employing Eqs. (|3"9")) -(RTjl). 

For convenience, here we use a rescaled concurrence 



a = (n- i)c, 



(51) 



and thus C r (0) = Cku(°) = Ct( )- Then we know that 
the two spin squeezing parameters and concurrence are 
same for the initial state. 



SPIN SQUEEZING AND CONCURRENCE 
UNDER DECOHERENCE 
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FIG. 1: (Color online) Time evolution of the spin squeezing 
parameters Cku> Ct> an d the rescaled concurrence CV for (a) 
N = 10 and (b) iV = 20. The inverse temperature is taken as 
P = 4/ojo. The insets in the figures show the magnification of 
the region where fjju nearly vanishes. The horizontal x axes 
are logarithmic, but the inset x axes are linear. 
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FIG. 2: (Color online) Time evolution of the rescaled concur- 
rence for different values of the inverse temperature p. Here, 
we choose TV = 10. 



The initial one-axis twisted state considered in this 
work is a symmetric state which can be expressed as a 
superposition of symmetric Dicke states. In other words, 
the N qubits behave effectively like a large spin N/2. 
After decoherence, not only the symmetric Dicke states 
will be populated, but also states with lower symme- 
try. Therefore, it is not sufficient to describe the sys- 
tem with only an (TV + l)-dimensional space. However, 



the exchange symmetry is not affected by the decoher- 
ence. In other words, a state with exchange symmetry 
does not necessarily belong to the maximally-symmetric 
space [l2; . Now by employing the hierarchy equation 
method, we calculate the spin squeezing parameters and 
the rescaled concurrence under decoherence, and com- 
pare the behaviors of them. 

As an example, we set the initial state given in Eq. (|4"T)l 
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50 100 150 200 

FIG. 3: (Color online) Time evolution of the spin squeezing 
parameter (| n for different values of /3, with TV = 10. The 
inset shows the magnification of the region when t is large. 

with 9 = 7r/10. The parameters of the Drude-Lorentz 
spectrum in Eq. pip are chosen to be A = 0.03ujq and 7 = 
0.15wo- In this section, we study the effect of the particle 
number TV and bath temperature T on the dynamics of 
spin squeezing and concurrence. 

Figures [Ul andQJ) show the time evolution of Cku> Ct 
and C r with two different particle number TV = 10 and 
TV = 20. The inverse temperature is set to be (3 — A/luq. 
The figures show that the decay rate of C r increases with 
TV. Although the rescaled concurrence of the initial state 
for TV = 20 is larger than that for TV = 10, it vanishes ear- 
lier. Also, the revival, after a sudden vanishing, becomes 
weaker with increasing TV. Both and Ct decay in an 
oscillatory way. We observe that £ T vanishes suddenly, 
while interestingly, £ku decays to zero asymptotically as 
shown in the insets. Comparing FigQJi andQJi , we find 
that for spin squeezing, the vanishing time changes little 
with increasing TV. 

Now we focus on the effects of the bath temperature 
on the dynamics of spin squeezing and rescaled concur- 
rence, which are shown by Figs. [2J3J These figures are 
plotted with a fixed particle number TV — 10 and dif- 
ferent temperature T . Here we choose the inverse tem- 
perature P = 4/wo, 3/wo, 2.5/wo) 2/wo) and we specially 
take P — 0.5/wo for t^u- Firstly, let us discuss the time 
evolutions of C r which are shown in Fig. [2j As expected, 
C r is suppressed with increasing temperature. When we 
choose a low temperature, such as /3 = 4/a>o, C r decays 
with multiple revivals. When the temperature increases, 
the revivals become weaker. C r even vanishes completely 
without revival when /? = 1/wrj- 

The spin squeezing is also suppressed with increas- 
ing T. As shown in Fig. [3j C^\j decays without sudden 
vanishing and approaches zero asymptotically (t — > 00) 
when temperature is not high enough, which is shown 
in the inset. Interestingly, when temperature reaches to 
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FIG. 4: (Color online) Time evolution of the spin squeezing 
parameter f° r different values of /3, with TV = 10. 



P = 0.5/ojo, Cku decays to zero quickly and suddenly 
without revival. The behavior is quite different with C r . 
While C T decays and suddenly vanishes even with low 
temperature as shown in Fig. 01 which is similar to C r . 

From the comparison, we find that the spin squeezing is 
not a satisfactory indicator of pairwise entanglement un- 
der decoherence for the open systemr, although they have 
close relations. Also, it is noted that the spin squeezing 
and concurrence both decay with oscillations, which is a 
reflection of the non-Markovian dynamics of the system. 



VI. CONCLUSION 

In this work, we consider an ensemble of N spin-1/2 
particles interacting with identical independent bosonic 
heat baths. The one-axis twisted state is chosen to be the 
initial state. The mean spin direction of the initial state 
is along the z-axis, and it does not change during the 
decoherence dynamics. For the open system we consider, 
we proved that the multi-qubit dynamics can be reduced 
into a two-qubit one. Then we use the hierarchy equation 
method to study the spin squeezing and concurrence un- 
der decoherence. This is an exact method without using 
rotating-wave and the Born-Markov approximation. 

From the numerical results, we find that the decay rate 
of the rescaled concurrence increases with the particle 
number N as well as the bath temperature T, and the 
revivals become weaker over time. For the spin squeezing, 
it is suppressed with increasing temperature as expected, 
while the vanishing time changes little with TV. The spin 
squeezing parameter (^yj vanishes asymptotically with 
low bath temperature and disappear suddenly when bath 
temperature is high enough. Interestingly, £ T vanishes 
suddenly even when bath temperature is low, which is 
similar to C r . 
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Appendix A: A formula on the partial trace 

It should be noted that the second term of Eq. (|2"91 is 

obtained by moving Yli=a u \ M horn the right of ptot (0) 
to the left, which use the relation of Eq. ff5U|) . We now 



7 

prove this property of partial trace. Considering a two- 
subspace case as an example: 

Tr 2 [{A l ®A 2 ) Pl2 {B 1 ®B 2 )] 
= Pi ® (("I A *) P12 (Si ® {B 2 |n)))] 

n 

= J2 i A i ® (("I Mp^Bx) ® (|m) (m| B 2 |n))] 

nm 

= i( A i ® (H B 2 |n) (n| A 2 ) pu-Bi ® |m)] 
= J2[A 1 ®((m\B 2 A 2 )p 12 B 1 ®\m)] 

m 

= r Tr 2 [A 1 ®{B 2 A 2Pl2 (B l ®I 2 )], (Al) 
where I\ 2 is the identity matrix of the 1 or 2 subspace. 
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